Rabi oscillations from ultrasound in spin systems 
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It is shown that ultrasound in the GHz range can generate space-time Rabi oscillations between 
spin states of molecular magnets. We compute dynamics of the magnetization generated by surface 
acoustic waves and discuss conditions under which this novel quantum effect can be observed. 
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Quantum mechanics of a spin cluster (e.g., a magnetic 
molecule) embedded in a solid is determined by the crys- 
tal field. The latter depends on the symmetry of the 
cluster and its environment Crystal field Hamilto- 
nians provide good description of molecular magnets at 
low temperature. For instance, dynamics of a spin S that 
prefers to look up or down along the anisotropy axis of the 
cluster can be described by a Hamiltonian H. = —DSl, 
where the anisotropy constant D arises from spin-orbit 
interactions. We are interested in the effect of the me- 
chanical rotation of the crystal on the molecular spin. 
For a similar problem involving an orbital moment, L, it 
is well known from classical mechanics that in the refer- 
ence frame rotating at an angular velocity il the Hamil- 
tonian acquires a term —KL ■ J7 (we use dimensionless L 
and S) . The same rule is expressed by the Larmor theo- 
rem in classical electrodynamics: Rotation is equivalent 
to the magnetic field, leading to the effective Zeeman 
term, — M • B, in the rotating-frame Hamiltonian, with 
M = ?i7L, B — 11/7, and 7 being the gyromagnetic 
ratio. The extension of the Larmor theorem to a spin 
is a consequence of the fact that in relativistic quantum 
theory the generator of rotations is J = L + S. Rigorous 
derivation of the term — (L-l-S)-O in the Hamiltonian can 
be obtained, e.g., from the study of the non-relativistic 
limit of the Dirac equation in the rotating frame 0] ■ 

Equivalence of the rotation to the magnetic field ex- 
plains Barnett effect Q : Rotation of a body of the mag- 
netic susceptibility x a-n angular velocity r2 generates a 
magnetic moment M — xO/7. The "spin- rotation cou- 
pling", — ftS • 17, can also lead to non-trivial quantum 
effects. Consider, e.g., a spin cluster with the Hamil- 
tonian Ti, ~ —DS^ that preserves the direction of the 
spin along the anisotropy axis Z due to commutation 
of Ti. with Sz- In the presence of the rotation about, 
e.g., the X-axis of the crystal the Hamiltonian in the 
rotating frame becomes Ti.' — ~DS^ — hSx^- This 
Hamiltonian, unlike Ti., does not commute with Sz and, 
therefore, allows transitions between the two orientations 
of S along the anisotropy axis. Thus, rotation alone 
can induce quantum transitions between spin states that 
are prohibited by the Hamiltonian of a stationary sys- 
tem H. We should emphasize that switching from the 



laboratory-frame Hamiltonian, Ti., to the rotating-frame 
Hamiltonian, Ti' — Ti — hS ■ CI, does not introduce any 
new spin-lattice interactions in addition to the crystal 
field. It is just another method to obtain solution of 
the problem, which, in the laboratory frame, requires in- 
troduction of the time dependence of the crystal field: 
e.g., Ti — —DSz, in the presence of rotation, becomes 
Ti — —D[n{t) ■ S]^ with n{t) being the instantaneous di- 
rection of the anisotropy axis. 

So far, quantum spin-rotation effects received little at- 
tention because only a very tiny magnetic field due to ro- 
tation can be produced in the rotating frame of a macro- 
scopic body. Consequently, the corresponding quantum 
effects have very low probability. This Letter is based 
upon the observation that local rotations of the crystal 
lattice produced by high-frequency ultrasound can eas- 
ily provide 10 G - 100 G fields in the rotating frame of a 
rigid spin cluster in a solid. Indeed, in the presence of 
the phonon displacement field, u(r, i), the angle of the 
local rotation of the crystal lattice, (50(r, t), and the cor- 
responding angular velocity, fl{r,t), are given by Q 



6cl){r,t) = ivxu(r,t), 



1 



n{r,t) = -Vxu(r,i) . (1) 



For a transverse sound wave of frequency / ~ 3 GHz and 
amplitude uq ^ Inm, this gives B — $1/7 ~ 10 G in 
the rotating frame coupled to the local crystallographic 
axes. Even greater local fields can be achieved with sur- 
face acoustic waves that have been recently used in ex- 
periments on molecular magnets ii. 

The equivalence of the effect of high-frequency trans- 
verse acoustic waves to the effect of high-amplitude ac 
magnetic field on paramagnetic spins immediately sug- 
gests that one can try to generate Rabi spin oscillations 
with the help of high-frequency ultrasound. Rabi ef- 
fect corresponds to the oscillation of the occupation 
numbers of two quantum levels in the presence of an ac 
field which frequency is close to the distance between 
the levels. On resonance, the frequency of Rabi oscil- 
lations is proportional to the amplitude of the ac field. 
The effort to observe Rabi oscillations between quantum 
states of molecular magnets in experiments employing ac 
magnetic fields [1, [^, liol . 11 1 has been going for some 
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time. For such experiments to succeed, the Rabi fre- 
quency must be greater than the spin decoherence rate. 
This typically requires the amplitude of the ac field to 
be greater than 1 G, which is not easy to achieve with 
electromagnetic waves but, as we have seen, is possible 
with surface acoustic waves. Note that the condition of 
the validity of the elastic theory, uq <^ X (where A is 
the phonon wavelength) automatically provides the con- 
dition <^ ll> = 27r/, which allows one to treat local 
rotations classically while treating the two-level system 
with level separation fko quantum-mechanically. 

For certainty we consider a crystal of molecular mag- 
nets with the anisotropy Haniiltonian 



Ha 



-DS^ + V, 



(2) 



where y is a small term responsible for the tunnel split- 
ting, A, of spin-up and spin-down states. The spin cluster 
is assumed to be more rigid than its elastic environment, 
so that the long-wave crystal deformations can only ro- 
tate it as a whole but cannot change its inner struc- 
ture responsible for the parameters of the Hamiltonian 
Ti.A- This approximation should apply to many molecu- 
lar magnets as they typically have a compact magnetic 
core inside a large unit cell of the crystal. We choose ge- 
ometry in which surface acoustic waves are running along 
the X-axis with the solid extending towards y > 0, see 
Fig. [TJ Using standard formulas ^] for the displacement 




FIG. 1: (color online) Geometry of the problem. 

field in a surface acoustic wave, one obtains 

5</)(r) = ^—uq e^*^'^ cos(kx — uit) Gz = 6(l>{x, t) , (3) 
2 ct 



where oj = Ctk£,, kt = k^\ — ^ is a real number 
between and 1 satisfying 



0, (4) 



and Ct^i are velocities of transverse and longitudinal 
sound. 

In the presence of deformations of the crystal lattice, 
local anisotropy axes defined by the crystal field are ro- 
tated by the angle given by Eq. ([3]). In the Hamiltonian, 



this rotation is equivalent to the rotation of the operator 
S in the opposite direction, which can be performed by 
the (25* + 1) X (25* + 1) matrix in the spin space 

S^i?-iSi?, R^e^-^'^. (5) 
The spin Hamiltonian in the laboratory frame becomes 



n 



(6) 



In order to find the laboratory- frame wave function j^*), 
it is useful to introduce the lattice-frame wave function 
l^(iat)^^ defined through the unitary transformation 



|^(/at)x 



*) 



(7) 



Differentiating it on time it is easy to see that this func- 
tion satisfies Schrodinger equation with the lattice-frame 
Hamiltonian 



where 



n 



d6(f) 



7:—uoe 
2 Ct 



smikx — tot) Gz 



(8) 



(9) 



To this point we have not made any assumptions about 
the magnitude ofScf), so that the equations ([6]) and ([8]) are 
exact. An interesting observation for the comparison of 
the effects of ultrasound and ac magnetic field is that the 
Hamiltonian ^ resembles the Hamiltonian of a particle 
of spin S in the ac magnetic field which amplitude scales 
as the square of the frequency. 

We are interested in the Rabi oscillations between the 
two lowest states of Ha- 



\^^) = —(\S)±\-S)), 



(10) 



where | ± S) satisfy S'^l ± 5) = iS*] ± 5). It makes 
sense, therefore, to project our Hamiltonian on the | ± S) 
states, making the problem essentially a two-state prob- 
lem. This gives 



■ (j\ — hujii sin(fca; — cut) 173 , 



(11) 



where A is the energy distance between the ground state 
10+) and the first excited state a^i = |S')(— 5| + | — 

S){S\, as = \S){S\ - \ - S){-S\, and 



2& 



(12) 



The two-state approach will be valid if A and hu are 
small in comparison with the distances to other spin lev- 
els. Note that the tunnel splitting A originates from the 
term V in Ha that does not commute with Sz- 



It is easy to check that a.t uj ^/fi, which is our case of 
interest for consideration of Rabi oscillations, the second 
term in Eq. (fTT|) can be treated as a perturbation as long 
as the wavelength of the acoustic wave satisfies A ^ Suq. 
For a not very large S this condition is always fulfilled 
by surface acoustic waves. The unperturbed eigenstates 
of the problem are then the eigenstates of cti given by 
Eq. (fTO|) . Their energies are ±A/2. The time-dependent 
perturbation produces transitions between these states, 
resulting in the Rabi oscillations when ~ A. The 
standard way to obtain the evolution of the wave function 
is to apply the rotating wave approximation 0]. Note 
that the coordinates x and y in Eq. can be viewed 
as parameters. Expressing the wave function as 



|*(i)("^*))=C+W|0+)+C_(O|^- 



(13) 



and starting with C_(0) = 0, C+(0) = 1 at t = 0,x = 0, 
one obtains 



C-{t) = ^e-*-sm 



. f^Rt 



(14) 



c+it) 

where 



fnRt\ A/h-uj 



(15) 



Assuming that every spin was in the ground state 10+) 
before the sound wave arrived, the spatial dependence 
of the wave function can be obtained by making a re- 
placement t ^ t — kx / (jj m. Eq. In the absence of 
spatial derivatives in the Hamiltonian (fTT|) . |\I'(t)^'''*)) is 
defined up to the phase factor exp[i0(a;, y)] with 9 be- 
ing an arbitrary real function of coordinates. From Eq. 
([7]) the wave function of the system in the laboratory 
frame is I*) = e~*'''^'^|^''^'°*)). Because of the smallness 
of 50, the dynamics of |^) essentially coincides with the 
dynamics of | 'I' ('"*)) and is given by Rabi oscillations be- 
tween the states \4>±) at the frequency Q.r. This is con- 
firmed by numerical calculations with the lattice-frame 
and laboratory- frame Hamiltonians, see Fig. [2] 

The expectation value of the projection of the spin onto 
the anisotropy axis (the Z-axis) is given by 



{m(t)\SA^{t))^2s'^x 



svt\{ijjt — kx) sin^ 



[Vint - Krx) 



+ — riji cos(cjt — kx) sin {Qnt — Krx) 



(16) 



where Kr = {VlR/Lo)k <C k can be called the "Rabi" 
wave vector. Thus, the space-time Rabi oscillations of 
the occupation numbers of spin states generate space- 
time oscillations of the magnetization of the crystal. On 




FIG. 2: (color online) Time dependence {t' = tA/h) of the 
probability to find the spin in the state a.t x = 0, S = 10, 
ujR = O.Ill), and u = 0.9(A//i). Dotted line (red): Numeri- 
cal result for the laboratory- frame Hamiltonian (|19p . Dash 
line (blue): Numerical result for the lattice- frame Hamilto- 
nian ITT]). Solid line (black): Analytical result given by Eq. 
(HI. 



resonance, when oj — A/h, Eq. (|16p simplifies to 

{^{t)\S^\'^{t)) ^ S cos{kx- tot) sin (ujRt-kRx) , (17) 

with kR = {u!R/uj)k <C k. The condition ujr <C w {Suq ^ 
A) implies that the time dependence of (Sz) at any point 
in space consists of the oscillations at frequency uj with 
beats of frequency cor. Similarly, (Sz) at any moment 
of time oscillates in space with the wave vector k and 
exhibits beats with the wave vector kR k. 

Our conclusions can be checked by obtaining the full 
solution of the problem in the laboratory frame in a par- 
ticular case of a biaxial symmetry, when V in Eq. (P]) 
equals E{Sl — Sy). In this case Eq. ^ reduces to 



H = -DSl + I |5ie-2»^^(-^*) + Si 



,2i5(t>{x.,t) 



}, (18) 



where S± = Sx i iSy. The second term can be treated 
as a perturbation provided that E ^ D. At cj <C 
(25+ l)D/h the dynamics of the wave function involves 
only a superposition of the \ iS) states. As in the lattice- 
frame consideration, it is then convenient to project the 
Hamiltonian p8)) onto these two states. In order to 
obtain such an effective two-state Hamiltonian that ac- 
counts for the tunnel splitting of the lowest energy states, 
one must apply perturbation theory for the degenerate 
states I ± S') to the S'-th order |l| . This results in 



-2iS<50(x,t) 



' ^{E/8D)'^ 



S){S\}. 
(19) 

is the tunnel- 



Here A = 8£'(2S')![(S'- 1)!] 
splitting for the biaxial model in the absence of lattice 
distortions [l^. Numerical solution for (Sz) that follows 
from Eq. (|19p . and its comparison with the analytical 
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Sz 




FIG. 3: (color online) Time dependence {t' = tA/h) of the ex- 
pectation value of the projection of the spin on the anisotropy 
a:xis at x — 0, S = 10, ujr — O.Ilj and uj — 0.9(A/ft). Dotted 
line (red): Numerical result for the laboratory-frame Hamil- 
tonian (|19ll . Dash line (blue): Numerical result for the lattice- 
frame Hamiltonian (lllf) . Solid line (black): Analytical result 
given by Eq. (fT6)l . 



higher frequencies of the acoustic waves will be required. 
In principle, surface acoustic waves of frequency as high 
as 100 GHz have been generated in experiment 1^. How- 
ever, since cor oc w^, raising to significantly may eventu- 
ally violate the condition <C lo under which our results 
were derived. 

In Conclusion, we have shown that transversal acoustic 
waves in the GHz range provide spin-rotation coupling 
that can be used to generate space-time Rabi oscillations 
in molecular magnets. When frequency of ultrasound, 
u), equals the distance between tunnel-split spin states, 
the magnetization on the surface of the crystal oscillates 
as (Sz) = S'cos(fca; — ujt) sm[uJiit — knx], where ojr = 
uPuqS/ {2ct) and kji = {ujj^/uj)k, with uq and Ct being 
the amplitude and the speed of the sound respectively. 

The authors thank Dmitry Garanin, Javier Tejada, 
Paulo Santos, and Jaroslav Albert for helpful discussions. 
This work has been supported by the NSF Grant No. 
EIA-0310517. 



solution given by Eq. are illustrated in Fig. [31 The 
beats discussed above are clearly seen in the figure. 

We shall now discuss conditions under which the above 
effects can be observed. The first condition is that the 
rate of decoherence of the spin states is lower than the 
frequencies involved. The lowest of these frequencies is 
lor = 2t:'^ PuqS / Ct, where /, uq, and Ct are the frequency, 
the amplitude, and the velocity of the sound. For, e.g., 
/ = 3 GHz, uo = 1 nm, S = 10, and ct = lO'^m/s, one ob- 
tains OJR'^ 2 GHz < w = 27r/ ~ 20 GHz. If such a high 
value of uor were to be produced by an electromagnetic 
wave, it would require the ac magnetic field of amplitude 
10 G, which is not easy to achieve in experiment. Note, 
however, that the Rabi oscillations oi {Sz) generated in a 
crystal of molecular magnets by ultrasound, contrary to 
the Rabi oscillations generated in a small crystal by an 
electromagnetic wave, will have a pronounced wave de- 
pendence on coordinates so that {Sz) averaged over the 
wavelength of the sound. A, will be zero. Consequently, 
measurements of the oscillations of {Sz) should be done 
on the scale that is small compared to A. 

Another restriction comes from the inevitable presence 
of the dc magnetic fields that generate the Zeeman en- 
ergy bias for the | ± S) states. Such fields can be of dipo- 
lar origin or they can be any stray fields in the system. 
They are not likely to aff'ect our results qualitatively if 
the Zeeman energy bias is small compared to ~ A. 
Note that the tunnel splitting. A, can be controlled by 
a transverse magnetic field. Thus, the above condition 
translates into Bi < 2TTf/{'jS) for the longitudinal field 
Bi. For / - 3 GHz one then needs Bi < 100 G. In the 
case of a greater bias field (and/or higher decoherence). 
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